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Abstract. The Planck mission will permit measurements of the polarization of the cosmic microwave background and of 
polarized foregrounds such as our own Galaxy with an unprecedented combination of accuracy and completeness. This will 
provide information on cosmological and galactic magnetic fields. The latter can be studied in detail via nearly Faraday- 
rotation free synchrotron and polarized dust emission. Methods are discussed to extract physically relevant information on 
the magnetic turbulence from Planck data and other measurements. 
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1. The Planck surveyor mission 
1.1. The experiment 

"Planck is a mission of the European Space Agency de- 
signed to answer key cosmological questions. Its ultimate 
goal is to determine the geometry and content of the Uni- 
verse, and which theories describing the birth and evolution 
of the Universe are correct. To achieve this ambitious objec- 
tive, it will observe the Cosmic Microwave Background ra- 
diation (CMB), emitted about 13 thousand million years ago, 
just 400,000 years after the Big Bang. Today the CMB perme- 
ates the Universe and is observed to have an average tempera- 
ture of 2.726 K. Small deviations from this average value (the 
so-called anisotropics), observable at angular scales larger 
than a few arcminutes, encode a wealth of information on 
the properties of the Universe in its infancy. The objective of 
Planck is to measure these properties with an unprecedented 
accuracy and level of detail." \The Planck Bluebook (2005)) . 

The Planck satellite will operate at the 2"'^ Lagrange point 
of the Sun-Earth system, with both Sun and Earth in the di- 
rection of the axis of the satellite, around which it rotates with 
1 rotation per minute. Two sets of detectors will be on board, 
the HEMT receivers for the Low Frequency Instrument (LFI) 
ranging from 30 to 70 GHz and bolometers for the High Fre- 
quency Instrument (HFI) ranging from 100 to 860 GHz. Both 
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instruments will have receivers sensitive to linear polariza- 
tion. 

The instrument beams will have resolutions of 5 to 30 
arcmin and observe the sky with an inclination of 85° with 
respect to the spin axis. About every 60 rotations, the spin 
axis will be rearranged to a new direction within a 10° cone 
pointing towards the Sun. Within seven months, the compli- 
cated scanning strategy of Planck will have covered the full 
sky at least once everywhere, and two complete coverages are 
planned. 



1.2. The data 

The CMB is not the only microwave emitter measured by the 
instruments. Galactic synchrotron, free-free, and dust emis- 
sion and extra-galactic sources contaminate the signal. For- 
tunately, nearly all of them have emission spectra which are 
very different from the CMB. The nine Planck spectral chan- 
nels cover one and a half orders of magnitude in frequency. 
This is necessary to separate the different physical compo- 
nents observed by a combination of spectral decomposition 
and spatial filtering. Detailed properties of the detectors can 
be found in Tabled The resulting physical component maps 
can then be analyzed separately according to their nature. 
Here, only the CMB and the Galactic emission components 
are of relevance. 



2 PRIMORDIAL MAGNETIC FIELDS 



2.2 Faraday rotation 



Table 1. Reproduced fmm \The Planck Bluebook (2005)| 

SUMMARY OF PLANCK INSTRUMENT CHARACTERISTICS 
LFI IIFI 

INSTRUMENT CHARACTERISTICS 



Detector Technology HEMT arrays Bolometer arrays 



Center Frequency [GHz] 


30 


44 


70 


100 


143 


217 


353 


545 


857 


Bandwidth (&v jv) 


0.2 


0.2 


0.2 


0.33 


0.33 


0.33 


0.33 


0.33 


0.33 


Angular Resolution (arcmin) 


33 


24 


14 


10 


7.1 


5.0 


5.0 


5.0 


5.0 


(5T/T per pixel (Stokes I)'' 


2.0 


2.7 


4.7 


2.5 


2.2 


4.8 


14.7 


147 


6700 


(5T/T per pixel (Stokes Q &U)'' 


2.8 


3.9 


6.7 


4.0 


4.2 


9.8 


29.8 







* Goal (in pK/K) for 14 months integration, 1 a, for square pixels whose sides are given in the row Angular Resolution. 



2. Primordial magnetic fields 

Particle physics scenarios exist which predict weak primor- 
dial magnetic fields, but usually on very small scales (see 
papers in these proceedings by M. Gasperini, by D. Dario, 
by K. Takahashi, and by D. Sokoloff). K. Subramanian and 
also T. Kahniashvili (these proceedings) take the pragmatic 
approach to assume that some primordial fields exist and to 
investigate its observable signatures. The same approach will 
be used in the following. 

2.1. Imprint in power spectra 

A uniform primordial magnetic field would break the cos- 
mic isotropy. Since this can be tested via CMB data, nG 
limits on a uniform field were derived from WMAP and 
other data, which certainly can be tightened by Planck (see 
[Giovannini (2005)l for a discussion). 

A tangled primordial magnetic field would lead to the 
presence of Alfven waves in the primordial plasma — pro- 
vided the field is dynamically significant and there is a non- 
empty scale interval between the scale of the field and the 
ion Larmor radius. The Alfvenic oscillations should imprint 
characteristic signatures onto the CMB. On large scales, the 
development of these oscillations is suppressed due to the 
very limited distance the relatively slow Alfven waves can 
have traveled for sub-equipartition fields (i? ^ 3 /^G from 
constraints of Big Bang nucleosynthesis). On small scales, 
"oscillations of ... Alfven waves get overdamped in the radi- 
ation diffusion regime, resulting in frozen-in magnetic field 
perturbations" ( |Jedamzik et al. (1998)^ . This freezing allows 
the magnetic fluctuations to survive until recombination and 
then to produce temperature fluctuations. 

[Subramanian & Barrow (1998)| estimate that fields with 
power-law spectra close to scale-invariance of 3 nG result 
in 10 /xK temperature fluctuations on multipole scales of 
I ~ 1000 - 3000 (below 10 arcmin). On sufficiently small 
scales, the magnetically induced power spectrum will easily 
be above the CMB spectrum of the acoustic oscillations due 
to Silk damping of the latter [Subramanian & Barrow (1998)| 
find that for the above parameters this happens at I ^ 2000. 

However, the identification of primordial fields using 
temperature fluctuations observed by Planck will be challeng- 
ing. Planck will have an angular resolution down to 5 arcmin 
(J. ^ 2000). Therefore, magnetic fields have to be stronger 



than assumed above in order for them to exceed the acous- 
tic fluctuations at these scales and be detectable. And even if 
Planck detects excess power at those scales, a large number 
of competing explanations is expected. 

A unique identification of CMB fluctuations due to mag- 
netic fields may be possible via the induced polarization sig- 
nature and/or non-Gaussian statistics. However, since these 
are more subtle measurements, their detection in Planck data 
is probably also very challenging. We (the authors) do not 
have a clear picture of the feasibility of uniquely identifying 
magnetic contributions to the CMB fluctuations. 

Nevertheless, Planck will be able to constrain primordial 
magnetic field scenarios since its measurements can always 
be used to set upper limits on magnetically induced fluctua- 
tions. 



2.2. Faraday rotation 

A nice signature of primordial magnetic fields would be an 
induced Faraday rotation of the CMB polarization. This ef- 
fect rotates polarization vectors with a characteristic de- 
pendence, which means that the lowest frequency channel 
of Planck will be best suited for such a signal (30 GHz). 
The technical specification of Planck requires an accuracy 
of polarization angle measurements of about 1 degree. This 
translates into a sensitivity in rotation measure of RM~ 
100rad/m2. 

The Faraday-rotation signal from nG magnetic fields is 
expected to be about one degree rotation at 30 GHz (see e.g., 
[Kosowsky et al. (2005)^ . However, since the CMB polariza- 
tion is a weak, small-scale signal, the changes due to Faraday 
rotation are even weaker and are also on small scales. The 
expected signal peaks with 0.1 /iK polarization fluctuations 
at I ~ 10''. These will be unobservable for Planck due to 
beam smearing effects (the 30 GHz system of Planck will 
be limited to 33 arcmin resolution I ^ 400), and due to the 
the overwhelmingly bright and highly polarized Galactic syn- 
chrotron emission at those frequencies extending to high lati- 
tudes. Usage of higher frequencies can decrease the Galactic 
synchrotron contamination at the price of an even stronger 
reduction of the wavelength-dependent Faraday -rotation sig- 
nal. 

In summary, a discovery of primordial magnetic fields via 
their Faraday effect is not to be expected from Planck. 



3.2 Reconstruction of large-scale fields 



3 GALACTIC MAGNETIC FIELDS 



3. Galactic magnetic fields 

The study of primordial magnetic fields with Planck will be 
difficult, partly because of the Galactic synchrotron contami- 
nation. However, the strength of the Galactic foreground sig- 
nal may open the possibility to study Galactic magnetic fields 
in great detail. 

3.1. Synchrotron, Faraday, and dust polarimetry 

Planck will be sensitive to magnetic fields B via 

1. synchrotron emission (B±), 

2. Faraday rotation (_B|[), 

3. and polarized dust emission {B±/Bj_), 

where ± and || refer to the field component with respect to 
the line of sight. 

Synchrotron emission is ideal for the study of galactic 
magnetic fields because the spatial distribution of relativis- 
tic electrons illuminating the fields should be smooth, due to 
their diffusivity, and can be measured by other means, e.g., 
via the inverse Compton scattering of CMB and starlight. 

Faraday rotation is also ideal because the required spa- 
tial distribution of thermal electrons can be measured by other 
means (e.g., their free-free emission). Despite the expected 
high precision in determining polarization angles (~ 1°), 
Planck will only be sensitive to RM > lOOrad/m^ due to 
its relatively high-frequency channels (> 30 GHz). Neverthe- 
less, Planck can provide "zero-wavelength" data for a multi- 
instrument Faraday campaign. 

Polarized dust emission due to scattered or partially ab- 
sorbed starlight (e.g.,jFosalb^^^(200^j^orspinning elon- 
gated dust grains (e.g. jDrain^^^azmarr(T998)| l might help 
to build and constrain models of the galactic field topology. 
However, it will be difficult to use this for quantitative in- 
vestigations of Galactic magnetic fields, due to the complex 
dependence of the signal on the poorly known dust properties 
and spatial distribution. 

3.2. Reconstruction of large-scale fields 

Given the large amount of information that Planck and other 
polarimetry instruments will provide on tracers of Galactic 
magnetic fields, one might ask if a full reconstruction of 
the large-scale field based purely on observational data is 
possible. The polarimetry of synchrotron emission provides 
three observables: the Stokes parameters /(A), Q(A), and 
C/(A), as functions of wavelength A. The goal would be to 
reconstruct from these the three magnetic-field components 
Bx{z), By{z), and Bz{z), which are functions of the 
coordinate along the line-of-sight z: 



The observed Stokes parameters are given in terms of 
their emissivities 
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dz £i{z, A), 



P{\) = Jdz ep{z, A) exp(2iA^ 0(z)), 

where P = Q + iU denotes the complex polarization, on 
which the rotation measure (f>{z) oc j dz Bz{z) ne{z) acts 
like a mathematical rotation operator. The emissivities are 
roughly given by combinations of the perpendicular magnetic 
field components: 

ei oc {BI + B^) TicRc, 

ep cc [Bl -Bl + 2tB,By) ncRe- 

The projection of the total intensity emission in the observa- 
tion removes any spatial information. For the polarized emis- 
sion, the frequency-dependent Faraday effect permits that 
spatial information imprints itself into the data. However, in- 
troducing the emissivity per Faraday depth 



ep{4>) = J dzep{z)S{(j){z) - (j)) 

permits the observable polarization to be expressed without 
any reference to the spatial distribution: 

/oo 
d4)£p{(t),\) exp(2a2(/>). 
-oo 

Therefore, a spatial reconstruction is impossible since spatial 
information is lost in the projection! 




Do we have enough information to reconstruct ep{4>)l 

P(A), expressed as a function of the variable a = 2 A^, 
and ep{(j)) are Fourier transforms of each other Knowing one 
function allows the reconstruction of the other. Unfortunately, 
measurements at negative values of a = 2 A^ would be re- 
quired to reconstruct the complex function ep{<f>), but this is 
impossible. Thus, at best half of the information necessary 
for a reconstruction of ep(0) is available. 

Although the above equations cannot be inverted, the 
forward approach is possible: the construction of polarized 
emission maps from Galactic models. 

The Fourier based formalism of Faraday rotation al- 
lows a computationally efficient calculation of spectral 
cubes of polarized maps using the fast Fourier transform. 
[WaeUcens (2005)| implemented the generation of maps of 
Stokes parameters I, Q, and U using a nested spherical pix- 
elization (HEALPix, |G6rski et al. (2005)| l in order to treat 
beam-depolarization effects accurately. This code requires as 
inputs 



Do we have enough information to reconstruct B{z)l 



- a Galactic electron model ( |Cordes & Lazio (2002)^ , 

- a Galactic cosmic -ray electron model, 

- a Galactic magnetic-field model. 



4 STUDYING MAGNETIC TURBULENCE 



4. 1 Faraday rotation & magnetic power spectra 




Fig. 1. Simulated Galactic total (top) and polarized (middle) 
intensity at 1.4 GHz and polarized intensity at the lowest 
Planck frequency of 30 GHz (bottom) calculated with the 
code of ]Waelkens (2005)] Note the Faraday depolarization ef- 
fects in the low frequency map that disappear at high frequen- 
cies. The units are erg/ (s cm^ Hz str). 

Fig- n shows synthetic emission maps produced by the code 
for a very simplistic magnetic field configuration (a logarith- 
mic spiral without field reversals but with an added random 
component). Surprisingly, the simulated maps already repro- 
duce many features of the observed ones, despite the simplic- 
ity of the model used. For instance, the observed Faraday- 
depolarization channels, which contribute to the low frac- 
tional polarization in the Galactic plane, seem also to be 
present in the simulation. 

4. Studying magnetic turbulence 

The Planck-polarization maps may allow a study of proper- 
ties of MHD turbulence in the Galaxy with a spatial resolu- 
tion not achievable in numerical simulations and a richness in 
detail not reachable in analytical investigations. 



But how to extract physically meaningful signals? In or- 
der to better understand the technical details of signal extrac- 
tion, we explain first, how magnetic power spectra could be 
extracted from Faraday rotation maps. Then we show that 
higher-order statistical properties of the magnetic turbulence 
may be measurable from Planck synchrotron-polarization 
data. 

4.1. Faraday rotation & magnetic power spectra 

Extended and polarized radio sources in galaxy clusters can 
be used to study magnetic turbulence in the intergalactic 
medium. The Faraday rotation map of Hydra A was analyzed 
by |Vogt & Enfilin (2005)| to determine the magnetic power 
spectrum within the Faraday screen in front of it (see Fig. 
O. This analysis was based on 

- homogeneous, isotropic, and divergence-free magnetic 
turbulence, 

- known cluster geometry (window function), and 

- Bayes' theorem connecting probabilistically the model 
spectra and the observed RM fluctuations. 

The incorporation of these assumptions in the analysis is ex- 
plained in the following. 

Homogeneous magnetic turbulence is best studied using 
the magnetic-field correlation tensor 

M,,(r) = {B,{x)Bj{x + r)), 

whose Fourier space representation is 

In general, this tensor would be described by nine functions 
defined in the three-dimensional fc-space. However, assuming 
isotropy and V .B = 0, it reduces to 



1 



— le 



H{k), 



2 V"'-' fc2 J k 
which depends only on two functions of fc = |fc|: 
3-d power spectrum: w[k) = y{B{k)-B{k)) 

3-d magnetic helicity: H{k) = j^^{B{k)-{B{k) x k)). 

Faraday rotation measures the line-of-sight projected 
magnetic field 



RM(a;^) = ao / dz ne{x) Bz{x) , 

J z^{x±) 

where ao = e'^/27rTO^c^. The RM autocorrelation function 

CRM(r±) = (RMix^) RM{x^ + r^)) 

is, therefore, connected to the magnetic-field correlation ten- 
sor Although we would need two measurable correlation 
functions in order to reconstruct the full tensor, the symmet- 
ric part of the tensor is fully encoded in the data (Enfilin & 
Vogt 2003): 



obs 
RM 



Here n^o is a characteristic electron density within the Fara- 
day screen of depth L. Since the information contained in the 
power spectrum alone (and not the helicity) determines the 
signal, the spectrum can be inferred from a single observable. 



4.2 Stokes correlators and Lorentz forces 



4 STUDYING MAGNETIC TURBULENCE 



For real data, the situation is a bit more complicated, since 
the limited size of the tested volume imprints on the mea- 
sured correlation function. This can be described by a win- 
dow function W{x), which contains the geometrical infor- 
mation on the variation of the cluster electron density and the 
limited lateral size of the radio source image. The connection 
between magnetic-field and RM spectra becomes: 

This integral equation translating w{q) into C^^{k±) has 
to be inverted in order to extract the former from the latter 
Bayes' theorem states that the probability for a model (here 
the power-spectrum 'w{q)) can be derived from the probabil- 
ity for the data (here the RM-map) given the model: 
P (model I data) oc P(data|modcl) P(model). 
Maximizing this probability with respect to the model param- 
eters provides the Maximum Likelihood power spectrum es- 
timator. This technique 

- is successfully used in CMB science and was tested for 
the Faraday application with mock RM data, 

- requires Gaussianity of RM fluctuations, but this is ex- 
pected due to the central limit theorem and is indeed in- 
dicated by observations, 

- takes care of the influence of the window function, and 

- provides errors and an error covariance matrix. 

It also provides a general conceptual blueprint for extracting 
statistics of physically meaningful quantities from the statis- 
tics of the observed signal. 

4.2. Stokes correlators and Lorentz forces 

Planck will provide high-accuracy and Faraday-free polariza- 
tion maps. Are physically meaningful quantities encoded in 
them? 

lEilek (1989)1 proposed to use correlation functions of 
Stokes parameters (henceforth Stokes correlators) to mea- 
sure magnetic power spectra. Since the power spectra are 
two-point second-order statistics, while the Stokes correla- 
tors are two-point fourth-order statistics, the latter cannot be 
translated into power spectra without some closure assump- 
tions. |Ejle^22l2l'^^'^°'^^'-'^^'-^'^ '■^^^ works for a Gaus- 
sian closure. However, MHD turbulence is not likely to be 
well represented by Gaussian fluctuations. This is evident 
from the appearance of coherent structures like flux sheets 
and filaments in numerical simulations (Fig.|3} and from the 
observation of filamentary magnetic fields in galaxy clusters 
(Clarke & EnBlin, these proceedings). 

If we were to look for physically meaningful quantities 

that are directly contained in the polarized emission data, the 

most obvious candidate is the Lorentz force 

1 1 

-J X B = V— + — B VB. 

c Stt in 

Since the Lorentz force is a quadratic quantity, its correlation 
function is a two-point fourth-order statistic. In the above ex- 
pression, the first term on the right-hand side is the magnetic 
pressure force, the second term is the magnetic tension force 

F^^B VB, 
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Fig. 2. Top: Faraday rotation map of the northern lobe of 
Hydra A (data: [Taylor et al. (I990)| map: |Vogt et al. (2005)| l. 
Bottom: Faraday-based magnetic power spectrum of Hydra 
A cool core ( |Vogt & EnBlin (2005)> . 

which not only determines the dynamical response of curved 
magnetic fields on the plasma, but can also be used to di- 
agnose the structure of the field: the tension-force correlator 
measures the gradients of the field along itself. Thus, {F{x) ■ 
F{x + r)) can be used as one of the quantitative measures of 
the folded structure evident in Fig. |3] where the field varies 
across itself on a much larger scale than along itself. Further 
discussion, as well as numerical measurements of the tension- 
force statistics can be found in | Sche kochihinetaLj2 004)| 
(see also the review by lSchekochihin & Cowley (2005)> . We 
propose to study the statistics of the tension force using po- 
larized radio observations. 

The correlation tensor of the tension-force or its Fourier 
counterpart 

where Cij^mn{k) is the Fourier-space representation of the 
two-point fourth-order magnetic correlation tensor 

C\j,mn{r) = {B.,{x)Bj{x)B„i{x + r)Bn{x + r)). 

A general isotropic fourth-rank tensor depends on 26 scalar 
functions Ci,,,26{r) of the distance r = |r|. Fortunately, ten- 
sor symmetries allow this set to be reduced to seven unknown 
scalar functions that fully determine the two-point fourth- 
order statistics of the field. In Fourier space, this reads 

+ C3{k)hkjk mkn 
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where Ci,,j{k) are all real. Is it possible to reconstruct these 
seven functions from observations? The accessible observ- 
ables are the Stokes parameters. These are related to the mag- 
netic fields in the emission region for Faraday-free, steep 
spectrum (a —1) radio emission:^ 



poo 

I{x^) = / dz[Bl{x) + Bl[x)\ 

J Zn 



Qix^) = / dz[Blix)-Blix)] 

J Zq 

U{x^) = / dz2B.,{x)B.y{x), 

J Zq 

from which we can construct six scalar Stokes correlators: 
in Fourier space, ^//(fej,), ilQQ{ki_), tuu{k±), S/Q(fe±), 
^iu{k±), and Y,Qij{k±). These can be expressed as linear 
combinations of the scalar functions Ci,,,5, while Cq and 
remain indeterminable. It turns out that the Stokes correlators 
are not independent and only four scalar correlation functions 
are, in fact, available from observations. Thus, we are short 
by one such function to fully reconstruct Ci,,,^. 

Despite this scarcity of observable information, the power 
spectrum of the tension force happens to be fully observable! 
It is completely expressed in terms of the Stokes correlators: 



2tt 



dip S„(fci)+3SQQ(fci) 



- \^T.QQ{k±) - T.uuil«±)j cos 4^3 + 4E/q(A;_l) cos2v3 

where if is the angle between k± and the x axis. This demon- 
strates that physically relevant information on MHD turbu- 
lence is encoded in the Stokes correlators and, therefore, can 
be obtained from Planck polarization data! The feasibility of 
this approach and the reliability of the underlying assump- 
tions clearly require further careful investigation, which is 
currently underway. 

5. Conclusions 

1 . The Planck surveyor mission is a high-precision experi- 
ment to study cosmology with the CMB. 

2. However, an unambiguous detection with Planck of pri- 
mordial magnetic fields — both from their imprint in the 
CMB power spectra and from their Faraday rotation — 
will be extremely challenging. But interesting constraints 
should be possible. 

3. Galactic magnetic fields can well be studied by Planck 
via synchrotron, Faraday, and dust polarimetry, allowing 
a model-based reconstruction of large-scale fields. 




' Here constant factors converting field strength to radio emissiv- 
ity have been suppressed for simplicity of the calculations. 



Fig. 3. Cross section of the field strength in the saturated state 
of a simulation of homogeneous isotropic MHD turbulence 
(run B in |Schekochihin et al. {2W^ . 

4. Studying magnetic turbulence with Planck is promising: a 
new technique of Stokes correlators may allow us to mea- 
sure the Lorentz-force power spectra in MHD turbulence, 
similar to the Faraday-based magnetic power spectra es- 
timates. 
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